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Scanning  Impedance  of  Proximity  Coupled 
Rectangular  Microstrip  Antenna  Arrays 


1.  INTRODUCTION 

Microstrip  antenna  technology  has  become  versatile  and  reliable,  and  features  such 
as  low  cost,  conformality,  and  reproducibility  make  the  printed  circuit  antenna  an  attractive 
candidate  for  many  antenna  array  applications.  A  primary  limitation  of  the  conventional 
microstrip  element,  however,  is  an  inherently  narrow  bandwidth.  Although  this  can  be  improved 
by  increasing  the  substrate  thickness,  feed  radiation  and  reactive  mismatch  make 
direct-coupled  microstrip  feed  lines  or  coaxial  probes  impractical  for  substrates  thicker  than 
approximately  0.05  dielectric  wavelengths. 

As  an  alternative,  proximity  coupling  can  be  used  to  overcome  the  substrate 
thickness  limitations  imposed  by  direct-coupled  feeds.  Figure  1  shows  three  different 
microstrip  elements  that  use  proximity  coupling.  Element  A  consists  of  two  stacked  patches; 
the  lower  patch  is  fed  by  a  coaxial  probe,  and  the  upper  patch  is  excited  by  the 
electromagnetic  field  of  the  lower  element.  Element  B  also  has  two  stacked  patches,  with  the 
lower  element  fed  by  a  coplanar  microstrip  line.  Element  C  is  a  single  patch  that  couples 
electromagnetically  to  a  microstrip  line  embedded  in  the  dielectric  beneath  it.  These  feed 
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A 


B. 


C. 


Figure  1.  Proximity-Coupled  Element  Geometries 

types  make  wider  bandwidths  possible  in  two  ways.  First,  the  overall  dielectric  thickness  can 
be  increased  by  adding  the  upper  layer  without  creating  additional  radiation  from  the  feeds  on 
the  lower  layer.  The  additional  volume  displaced  by  the  element  translates  into  a 
lower  Q,  and  as  a  consequence,  the  bandwidth  is  increased.  Secondly,  for  elements  A  and  B, 
wider  bandwidths  are  obtained  by  offsetting  the  resonant  frequencies  of  the  two  patches. 
Experimental  results  for  isolated  proximity-coupled  microstrip  antenna  elements  with  bandwidths 
of  up  to  26  percent  have  been  reported  in  References.1’2,3’3 4  Although  these  results  indicate  a 
substantial  improvement  over  conventional  microstrip  antenna  bandwidths  of  5-10  percent,  no 
studies  have  been  done  to  assess  the  effects  of  mutual  coupling  and  dielectric  surface  waves 
in  a  large  array  environment. 

This  report  describes  a  full  wave  analysis  that  uses  the  infinite  array  spectral 


1.  Yasuo,  S.,  Miyano,  N.,  and  Chiba,  T.  (1983)  Expanding  the  Bandwidth  of  a  Microstrip  Antenna, 
IEEE  AP-S  Int.  Symp.  Digest,  pp  366-369. 

2.  Sabban,  A.  (1983)  A  New  Broadband  Stacked  Two-Layer  Microstrip  Antenna,  IEEE  APS  Int. 
Symp.  Digest,  pp.  63-66. 

3.  Chen,  H.,  Tulintseff,  A.,  and  Sorbello,  R.  (1984)  A  Broadband  Two-Layer  Microstrip  Antenna. 
IEEE  APS  Int.  Symp.  Digest,  pp.  251-254. 

4.  Cock,  R.  and  Christandoulous,  C.  (1987)  Design  of  a  Two-Layer,  Capacitively  Coupled, 

Microstrip  Patch  Antenna  Element  for  Broadband  Applications,  IEEE  APS  Int. 

Symp.  Digest,  pp.  936-939. 
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dyadic  Green’s  function  for  electric  current  sources  embedded  in  a  two-layer  grounded 
dielectric  slab.  This  approach  has  been  successfully  applied  to  several  printed  antenna 
geometries5,6, 7  and  has  the  advantage  that  it  includes  the  effects  of  surface  waves  and  replicates 
the  mutual  coupling  environment  seen  by  the  elements  in  a  large  phased  array. 

A  method  of  moments  is  used  to  solve  the  electric  field  integral  equation,  and 
three  types  of  expansion  modes  are  used  to  efficiently  model  the  current  distribution  on  the 
feed  line  and  patches.  A  technique  is  developed  for  modeling  microstrip  feed  lines  in  an 
infinite  array  unit  cell. 

Theoretical  results  are  presented  for  input  impedance  and  scanning  bandwidth  as 
functions  of  substrate  parameters,  array  spacing,  element  geometries,  and  feed  types.  The  scan 
performance  of  each  of  the  three  element  types  is  studied.  It  is  found  that  scan  blindness 
conditions  do  occur,  and  that  the  blindness  angle  moves  towards  broadside  with  increased 
dielectric  thickness  or  element  spacing.  A  complete  reflection  is  predicted  at  the  E-plane 
blindness  of  element  A,  and  a  large,  non-unity  reflection  is  predicted  for  elements  B  and  C. 
Experimental  data  from  a  waveguide  simulator  and  an  1 1  x  1 1  array  are  shown  to  agree  well  with 
the  results  computed  by  the  model. 

The  fundamental  bandwidth  limitations  of  proximity  coupled  array  elements  are 
investigated,  and  it  is  found  that  the  maximum  broadside  bandwidth  that  can  be  attained  by  the 
stacked  patch  element  is  approximately  25  percent.  This  limit  is  showi  to  be  related  to  the 
theoretical  limit  for  broadband  matching  of  a  complex  load  impedance,  as  determined  by 
Fano.5 6 7 8 

A  sample  array  element  is  designed  for  maximum  scanning  bandwidth,  and  it  is  shown 
that  there  is  a  tradeoff  between  bandwidth  and  scan  range.  The  sample  element  is  capable  of  a 
10  percent  bandwidth  (2:1  VSWR)  over  a  50  degree  scan  cone. 

The  computed  and  measured  scanning  impedance  results  indicate  that  mutual  coupling 
can  have  a  substantial  effect  upon  resonant  frequencies  and  bandwidth  in  a  scanning  array 
environment. 


5.  Pozar,  D.  and  Schaubert,  D.  (1984)  Scan  Blindness  in  Infinite  Arrays  of  Printed  Dipoles, 

IEEE  Trans.  Antennas  and  Prop.,  AP-32:602-610. 

6.  Pozar,  D.  and  Schaubert,  D.  (1984)  Analysis  of  an  Infinite  Array  of  Rectangular  Microstrip 
Patches  with  Idealized  Probe  Feeds,  IEEE  Trans.  Antennas  and  Prop.,  AP-32:1 101-1 107. 

7.  Pozar,  D.  (1989)  Analysis  of  an  Infinite  Phased  Array  of  Aperture  Coupled  Microstrip  Patches, 
IEEE  Trans.  Antennas  and  Prop.,  AP-37:41 8-425. 

8.  Fano,  R„  (1958)  Theoretical  Limitations  on  the  Broadband  Matching  of  Arbitrary  Impedances, 
Journal  of  the  Franklin  Institute,  249:57-84  and  139-154. 
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2.  THEORETICAL  ANALYSIS 
2.1.  Integral  Equation  Formulation 


This  work  uses  the  electronic  field  integration  formulation  to  study  the  input  impedance 
and  radiation  properties  of  proximity-coupled  microstrip  antenna  array  elements.  To  derive 
the  electric  field  integral  equation,  we  begin  with  Maxwell’s  equations  for  time  harmonic 
fields; 

V  x  E(x,y,z)  =  -  j(OJJ.H(x,y,z)  (1) 

V  x  H(x,y,z)  =  jcoeE(x,y,z)  +  J(x,y,z).  (2) 


An  expression  for  E  is  obtained  by  taking  the  curl  of  Eq.  (1)  and  substituting 
the  curl  of  H  on  the  right-hand  side  with  Eq.  (2); 

V  x  V  x  E(x,y,z)  -  k2E(x,y,z)  =  -jwpj(x,y,z)  (3) 

where  k  is  the  propagation  constant  defined  as; 

k2  =  <o2|i£.  (4) 

Eq.  (3)  relates  the  electric  field  E  directly  to  the  electric  current  density  J. 

Next,  define  a  dyadic  Green’s  function  to  be  the  solution  to  a  corresponding 
equation; 


V  x  V  x  G(x,y,z|x',y',z')  -  k2G(x,y,z| x'^'.z')  =  -I  5(x-x')8(y-y')5(z-z')  (5) 

where  I  is  the  identity  dyad  and  5  is  the  Dirac  delta  function.  Physically,  the  Green’s 
function  can  be  interpreted  as  the  electric  field  of  an  infinitesimal  electric  dipole  located 

at  x=x’,  y=y’,  z=z'.  Tn  dyadic  notation,  G  is  defined  as; 


AAAAAAAAAAAA 

G  =  xG  x  +  xG  y  +  xG  z  +  yG  x  +  yG  y  +  yG  z 
xx  xy3  xz  3  yx  3  yy3  3  yz 

A  A  A  AAA 

+  zG  x  +  zG  y  +  zG  z 
zx  zy3  zz 


4 


where  G  represents  the  i -directed  electric  field  of  a  i -directed 
>  j 

infinitesimal  electric 

dipole.  In  general,  all  "ine  components  of  the  dyad  are  required  to  completely  describe  the 
electric  field. 

The  electric  field  integral  equation  is  oh’dned  by  combining  Eq.  (3)  and  (5), 
applying  Green’s  second  identity,  and  choosing  appropriate  boundary  conditions  for  G  to  get; 


E(x,y,z)  =  jeofj. 


G(x',y',z')  •  J(x,y,z  I  x'.y'.zOdx'dy'dz'. 


(7) 


2.2.  Fourier  Transform  Green’s  Function 

The  integral  equation  formulation  for  the  analysis  of  proximity-coupled  microstrip 
array  elements  of  Figure  1  requires  the  Green’s  function  for  a  two-layer  grounded  dielectric 
slab.  As  mentioned  earlier,  this  is  identical  to  solving  for  the  fields  of  a  point  source 
embedded  in  two  layers  of  dielectric  media  above  a  conducting  plane,  as  shown  in  Figure  2, 
where  J  is  given  by; 


Figure  2.  Delta  Function  Current  Source  in  a  Two-Layer  Grounded  Dielectric  Slab 
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(B) 


J  =  x  5(x-x')5(y-y')5(z-h2) 

The  desired  fields  can  be  found  by  working  directly  with  the  wave  equations  for  E 

Q  ,  ^ 

and  H  ,  and  finding  the  transverse  fields  from  these  components.  The  wave  equations  in  a 


source-free  region  are; 

V  x  V  x  E  -  k2E  0 

(9) 

V  x  V  x  H  -  k2H  =  0 

(10) 

where  k2  =  ,  and  e  =  e2eo  for  0<z<h2,  e  =  for  h2<z<ht+h2  ,  and  e  =  eQ  for  z>h]+h2  . 

For  the  geometry  under  consideration,  the  boundary  conditions  occur  at  constant 
z-planes.  The  Fourier  transform  method  is  well  suited  for  solving  differential  equations  with 
boundary  conditions  that  are  unbounded  in  one  or  more  dimensions.  A  Fourier  transform  pair  is 
defined  as; 


E(x,y,z) 


oo  oo 


.00  -00 


(11) 


E(kx,ky,z) 


E(x,y,z)  e’-*^xx 


e-jy 


dx  dy. 


(12) 


.OO  .OO 

The  wave  Eqs.  (9)  and  (10)  are  then  Fourier  transformed  in  each  region,  thereby 
reducing  the  problem  to  one  spatial  dimension.  General  solutions  for  E^  and  H?  are  therefore; 


E  =  A 

-jk  z 
eJ  0 

z  >  h  +h 

(13) 

z 

1  2 

/ 

HZ=B 

-jk  z 
e  1  0 

z  >  h  +h 

1  2 

(14) 

9.  Pozar,  D.  (1989)  Analysis  and  Design  Considerations  for  Printed  Phased  Array  Antennas, 
Microstrip  Antenna  Handbook,  J.  James  and  P.  Hall,  Eds.,  Peter  Peregrinus,  London. 


6 


E_  =  C  e‘jk,z  +  D  e>k,2 


h~  <  z  <  h  +h 
2  1  2 


H  =  E  e‘jkiZ  +  F  e’kiz 


h  <  z  <  h  +h 
2  1  2 


E  =  G  e'jk2Z  +  H  e^k2Z  0  <  z  <  h 


H  =  I  e'Jk2Z  +  J  e>k2Z  0  <  z  <  h 

Z  2 


where 


k  2  =  k2  -  P2 
o  o  r 


k,2  -  e,k^  -  p2 


k  2  =  e  k2  -  p2 
2  2  0  r 


P2  =  k2  +  k2 
*  y 


In  the  transform  domain,  the  transverse  field  components  can  be  expressed  in  terms  of  E 


and  H  in  each  region9; 


E  -  a—  f  w4"ky  H 

X  ^  az  Z  p2  z 


jkv  a  «  wM-nky  ~ 

E  = - £  %-  E  +  — —  H  . 

y  p2  dz  Z  p2  z 


we  e  k  ~ 


Ez 


j  k  „  we  £  k  ~ 

u  _  y  o  fj  0  r  X  p 

Hy  ~ p2  dz  z  p2  Ez 
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These  equations  apply  to  all  three  regions,  with  set  accordingly  to  e  ,  e  ,  or  1.  The 
boundary  conditions  of  the  problem  are; 

1.  E  and  E  are  zero  at  z=0. 

*  y 

2.  E  and  E  are  continuous  at  z=h  . 

x  y  2 

3.  E  and  E  are  continuous  at  z=h  +h  . 

x  y  12 

4.  H  is  continuous  at  z=h  . 

x  2 

5.  H(lt(z=h  )  -  Hm(z=h  )  =  -e'Kk.xVk  y  ) 

6.  H  and  H  are  continuous  at  z=h  +h  . 

X  y  2  1 

Note  that  a  change  in  the  direction  of  the  infinitesimal  source  or  the  z-coordinate  of  the 

source  will  require  a  modification  to  conditions  4,  5,  and  6.  Eqs.  (12)  -  (17)  for  E  and 

H?  are  substituted  into  Eqs.  (21)  -  (24),  and  by  applying  the  boundary  conditions  above,  a  set  of 
10  equations  can  be  solved  algebraically  for  the  10  unknowns.  This  process  is  repeated  for 
each  of  the  three  source  components  x,  y,  and  z,  so  that  a  complete  field  dyad  is  obtained  for 
a  source  at  z=h2-  The  process  is  repeated  for  a  source  at  z=hf+h2.  We  shall  designate  the 

transformed  field  dyads  as  Q(k^,k^,z=h2)  and  Q(k^,k^,z=h]+h2),  and  they  represent  a  complete 
description  of  the  transformed  fields  of  arbitrarily  oriented  infinitesimal  dipoles  at  z=h2 

and  z=h(+h2.  Algebraic  expressions  for  these  Q’s  are  given  in  the  Appendix. 

To  recover  the  x  and  y  dependencies,  an  inverse  Fourier  transform  is  taken; 

oo  oo 

G(x,y  I  x',y')  =  j  j  ^.ky^^^'^Vy'^d^dky  (26) 

.OO  -OO 


where  G(x,ylx',y')  is  the  dyadic  Green’s  function  for  a  source  at  x=x’,  y=y’. 

The  spectral  dyadic  Green’s  function  of  Eq.  (26)  can  be  generalized  to  an  infinite 
array  of  sources.3  Figure  3  shows  the  geometry  of  an  infinite  array  of  horizontal  electric 
dipoles  with  periodic  spacing,  a  in  the  x-direction  and  b  in  the  y-direction.  The  p.q01  source 
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is  located  at; 


x  =  x'  +  pa 

p 

(27) 

yq  =  y'  +  qb. 

(28) 

z 


Figure  3.  Infinite  Array  of  Delta  Function  Current  Sources  in  a  Two- Layer  Grounded 
Dielectric  Slab 

In  order  to  scan  the  array  to  an  angle  (9,<{>),  the  currents  on  the  p.q*  source  must 
be  phased  as; 

e«w  _  e~iko(pau  +clbv)  (29) 

where 

u  =  sinG  cost))  (30) 

v  =  sinG  sin<|)  (31) 

The  dyadic  Green’s  function  of  the  infinite  array  is  then  given  by  the  doubly  infinite 
superposition  of  fields  from  each  infinitesimal  element; 
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oo 


©o 


(32) 


Substituting  in  the  expression  for  from  Eq.  (26); 


OO  oo 

x  y  x  y 

.00  .00 


(33) 


As  shown  in  reference  5,  the  Poisson  sum  formula  can  be  applied  to  Eq.  (33)  to  reduce  the 
doubly  infinite  summation  of  doubly  infinite  integrals  to  a  doubly  infinite  summation; 


00  OO 


0,(x.ylx'y,  =  j>  [  [  wx.k  Jk,(*-*vyy-y'> 

p  =  -oo  q  -  .00  * 

(34) 

where 

kx  =  JiE-+ko" 

(35) 

ky  =  -T9-  +  k»V- 

(36) 

Eqs.  (26)  and  (34)  for  the  single  and  infinite  array  Green’s  functions  can  then  be  used 
in  the  integral  equation  for  the  field  of  an  arbitrary  current  distribution  as  given  by  Eq.  (7) 
For  the  single  element,  Eq.  (7)  becomes; 


E(x,y,z) 


J(x',y\z'). 


00  00 

Q(k  ,k  ,z  I  z')e^x^x"x  ^e^y^'y  ^  dk  dx'dy'dz' 
x  y  x  y 

* 

.00  .00 


For  the  infinite  array,  Eq.  (7)  can  be  written  as; 


(37) 
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}(0\X 

E(x.y,z)  = 


J(x',y',z'). 

oo  oo 

I  I  Q(kx.k¥)eikjI<X'X,)e’ky<y‘yVd/d^  (38) 

P  =  -OO  q  =  -OO  * 

with  k  and  k  given  by  Eqs.  (35)  and  (36).  Note  that  the  current  distribution  J  in  Eq.  (38) 

*  y 

need  only  be  integrated  over  the  unit  cell  of  the  infinite  array. 

2.3  Moment  Method  Solution 

The  unknown  electric  current  density  J  can  be  solved  approximately  by  applying  the 
moment  method  procedure  to  the  electric  field  integral  in  Eqs.  (37)  and  (38).  The  boundary 
condition  that  the  total  tangential  electric  field  must  equal  zero  on  all  conducting  surfaces 
can  be  written  as; 

E  =  E™  +  E*4'  =  0.  (39) 

tan  tan  tan 

Using  Eq.  (7),  this  can  be  rewritten  as; 


J"0  •  G  dsinc  = 

* 

, 

, 

J  4C4‘-  G  ds 


scat 


S  i  nc 


S  scat 


(40) 


where  Jmc  is  the  incident  surface  current  density,  and  J**1  is  the  unknown  scattered  current 

density.  As  will  be  discussed  later,  Jmc  corresponds  to  the  impressed  current  density  on  the 
vertical  probe  for  element  A,  and  is  the  incident  traveling  wave  current  for  elements  B  and  C. 

J5""  is  the  surface  current  density  that  is  induced  on  the  feed  line  and  patches  by  the 
fields  of  the  incident  current  density. 

The  unknown  function  J*6*1  is  then  expanded  in  a  set  of  linearly  independent  basis 
functions  J  with  unknown  coefficients  I  ; 

n  n 

N 

J(x,y,z)  =  £  InTn(x,y,z).  (41) 

n  =  I 

The  expansion  modes  are  chosen  to  closely  resemble  the  anticipated  current  distributions. 
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Figure  4  shows  the  three  types  of  modes  that  are  used  to  analyze  the  proximity  coupled 
elements,  and  they  will  be  described  in  detail  in  the  next  section. 


Z 


✓  ENTIRE  DOMAIN 

TRAV.  WAVE 

~  "TVjVx5*?'*  < 

PWS 

€ 

2 

XB  XL 

Figure  4.  Expansion  Modes  for  Proximity-Coupled  Element  Geometries 


Eq.  (41)  is  then  substituted  into  Eq.  (40)  and  multiplied  by  a  testing 

function  J  and  both  sides  of  Eq.  (40)  are  integrated  over  the  conducting  surface  of  the 
testing  function,  s  ; 


Jmcds  ds  = 

m  inc 


J  •  G  •  J  ds  ds 

m  n  m  n 


S.  S 
me  m 


S  S 
n  m 


m  =  1,2, ...,N  ;  n  =  1,2,. ..,N 


(42) 


This  testing  procedure  leads  to  a  matrix  equation  that  can  be  solved  for  the  unknown  I  ’s; 


[Z 1  ‘[V]  =  [I] 


(43) 


where 


Z 


mn 


J  ds  ds 

n  m  n 


S  S 
n  m 


(44) 
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V  =  - 
m 


s .  s 

inc  m 


J  •  G  •  Jmc  ds  dsmc  . 

m  m 


For  the  single  element  case,  the  Green’s  function  of  Eq.  (26)  can  be  substituted  into  Eqs. 
(44)  and  (45)  to  get; 


,s ingle  _ 


J(x.y)  •  J  (x'y) 


f  CKki  dk  dx'dy'axdy 

X  y  x  y 


*  A 


J  (x>y)  •  J.  (x',y') 

m  me 


S  .  S 
j  nc  m 


Q(k  ,kv)ejkx(x'x,)ejky(y'y/)dk  dk  dx'dy'dxdy. 
x  y  x  y 


Next,  a  Fourier  transform  of  the  electric  surface  current  densities  can  be  defined  as; 


W9  '  JJ 


y)  e  ^kxX  e ’^kyy  dxdy. 


The  exponential  factors  of  the  transformed  Green’s  function  in  Eqs.  (46)  and  (47)  can  be 
incorporated  with  the  surface  integrals  to  form  Fourier  transforms  of  the  electric  surface 
current  densities,  and  the  expressions  for  Z  and  V  can  be  rewritten  as; 


mn  m 


ry%  ingle  _  1 

- 


i  r 

?  J  J 


F  (-k  .-k  )  •  Q(k  Jt )  ■  F  (k  Jc )  dk  dk 

mxy  xy  nxy  xy 
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The  Fourier  transform  integrations  in  Eq.  (48)  are  done  in  closed  form,  so  the  computation  of 
the  impedance  matrix  [Z]  and  voltage  vector  [V]  is  very  efficient. 

Similarly,  when  the  infinite  array  Green’s  function  of  Eq.  (34)  is  substituted  into 
Eqs.  (44)  and  (45),  the  expressions  for  Z  and  V  can  be  reduced  to; 

mn  m 


Vinf  =  '4  Y  Y  F  (kic)  •  Q(k  k  )  •  F.  (-k.-k) 

m  ab  L  L  m  x  y  ^  x  y  mcv  x  y 


p  =  -oo  q-.oo 


y  me  x  y 


with  k  and  k  given  by  Eqs.  (35)  and  (36). 

x  y 


(51) 

(52) 


2.4.  Expansion  and  Test  Modes 

On  the  rectangular  microstrip  patches,  entire  domain  basis  functions  are  used  because  of 

A 

their  correspondence  with  the  cavity-like  behavior  of  the  patches.  For  x-directed  currents, 
the  modal  functions  are  given  by; 

T^ch(x,y)  =  x  -L  sin[-~^-[x  -  ^-]]  (53) 

where  p  and  p  are  the  patch  dimensions,  and  r  is  an  integer  index  accounting  for  the  number 
*  y  a 

of  x-variations.  The  expression  for  the  y-directed  entire  domain  current  s  is  given  by; 

T^x.y)  =  y  -L  sin[-^-|y  -  (54) 

^  y 

where  s  is  an  integer  accounting  for  the  number  of  y- variations.  Note  that  the  functions  in 

A  A 

Eqs.  (53)  and  (54)  are  constant  in  the  y  and  x  directions,  respectively.  The  Fourier 

A 

transform  of  the  x-directed  modes  is  given  by; 

F^tch(kx,ky)  =  x  e'jkx(px/2)  e‘-’ky(Py/2) 
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(55) 


For  the  y-directed  modes,  the  Fourier  transform  is  given  by. 


F^k  ,k  )  =  y  e"jkx(p*/2)  e'jky(p/2) 

n  x  y 


[stt/p  ]  pkyPy  -  COs(sJt)| 

[  KT7^ 

k2 

X 

The  currents  on  the  microstrip  feed  lines  are  modeled  by  incident  and  reflected 
semi-infinite  traveling  wave  modes  and  sub-domain  piecewise  sinusoidal  modes  as  in  references 
10  and  11.  The  traveling  wave  modes  are  the  most  efficient  way  to  represent  the  quasi-TEM 
modes  that  propagate  on  a  microstrip  line,  and  the  piecewise  sinusoidal  modes  model  the 
non-propagating  currents  excited  by  the  discontinuity  in  the  microstrip  lines  of  elements  B 
and  C.  Overlapping  several  piecewise  sinusoidal  modes  on  the  lower  patch  of  element  B  helps  to 
enforce  continuity  of  current  at  the  line-patch  intersection.  More  details  on  the  feed 
modeling  will  be  given  in  the  next  section. 

The  incident  currents  are  given  by  the  expression; 


0*,y)  = 


-jk  x 


x  e 

a  -ik* 

x  e  cos 


fe>] 


<  x  <  X„ 


XD  <  x  <  X 

B  L 


,  -  W 
y  I  <  T 


I  y  I 


(57) 


where  k  is  the  effective  propagation  constant  of  a  quasi-TEM  mode  on  the  microstrip  line.  As 
seen  in  Figure  4,  the  traveling  wave  modes  terminate  with  a  one-quarter  cycle  of  a  cosine 
function  beginning  at  x=XQ  and  ending  at  x=XL>  This  feature  is  helpful  in  reproducing  the 
shape  of  the  current  distribution  at  the  end  of  the  microstrip  line,  and  was  found  to  reduce 
the  number  of  piecewise  sinusoidal  modes  needed  to  accurately  model  the  currents  there. 

The  transverse  current  distribution  is  assumed  to  be  constant,  and  this 
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approximation  has  been  found  to  be  valid  in  earlier  studies.10’11  A  proportionality  constant  of  1AV 
is  used  to  normalize  the  current  amplitude,  where  W  is  the  width  of  the  microstrip  line.  Using 
a  technique  developed  in  reference  10,  k^  is  found  by  using  the  electric  field  integral  Eq. 

(7)  for  the  field  of  an  infinite  microstrip  line,  and  enforcing  the  boundary  condition  that 
the  total  tangential  electric  field  is  equal  to  zero  across  the  microstrip  line.  By  applying 
the  Fourier  transform  method  described  earlier  with  a  traveling  wave  mode  on  the  microstrip 
line,  k^can  be  computed  as  the  root  of  the  characteristic  equation; 

oo 

K  V  dky  = 0  (58) 

-OO 


where 

F  (k)  = 
y  y 

is  the  Fourier  transform  of  the  transverse  current  distribution,  and  Q  is  given  in  the 

XX 1 

Appendix. 

The  reflected  traveling  wave  currents  are  given  by; 


(59) 


10.  Jackson,  R.  and  Pozar,  D.  (1985)  Full- Wave  Analysis  of  Microstrip  Open-End  and  Gap 
Discontinuities,  IEEE  Trans.  Microwave  Theory  and  Tech.,  MTT-33: 1036- 1042. 

11.  Voda,  S.  and  Pozar,  D.  (1987)  A  Rigorous  Analysis  of  a  Microstrip  Line  Fed  Patch  Antenna, 
IEEE  Trans.  Antennas  and  Prop.,  AP-35: 1343- 1350. 
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W*?]  -~j(W 

[2  <  x"-x8  >  ]  -  (k.+kJ 

The  Fourier  transform  of  the  reflected  travelling  wave  mode  is  given  by; 


F™(k  ,k)  =  x  F  (k  )  ■  TtSfk-kl 

ref  X  y  y  y  [x  ej 


e'^'V^ 


je-j<k,-ke)XB 


k  -k 


■)X|-[irx^] +  j(kA) 
[’(x"  x — s]  -  (k,-  kJ 


-j(k  -k  )x 
e JV  x  t  b 


The  piecewise  sinusoidal  modes  are  expressed  mathematically  as; 


TPWS,  x  A 

\  (x,y)  =  x 


si"[krr~“  lx_xJ]  U-xJ  <  4 

w  sin  ^ke-4y— |  I  y  I  <  -  y-- 


ly  I  < 


where  L  is  the  length  of  the  mode,  and  Xn  is  the  x-coordinate  at  the  center  of  the  mode.  The 
Fourier  transforms  of  the  piecewise  sinusoidal  modes  are  given  by; 


■  j  *  *  r  i  1 

_pWS  A  e  cos  k  -y-  -  cos  k 

F^k  ,k  )  =  x - — - 1/  2  J - p 

n  x  y  sin  k  k2  -  k2 

e  L  Ox 


F  (k  ) 
y  y 


2.5.  Scanning  Impedance  Analysis 

The  modeling  techniques  described  in  the  previous  sections  can  be  customized  to 
study  each  of  the  three  proximity  coupled  element  geometries  in  Figure  1.  The  vertical  probe 
feed  of  element  A  is  modeled  by  an  infinitesimal  electric  current  filament  at  x=x  ,  y=y  ; 

p  p 

T:;>y)  =  *  5<x-xp)5(y-yp)  (65) 

The  Fourier  transform  of  this  function  is  simply; 
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(66) 


F"^  ,k  )  =  i  e-Sk„\  e'*/, 

inc  x  y 

The  voltage  vector  elements  are  computed  by  means  of  the  discrete  spectral  series  summations 
of  Eq.  (52); 


V  "f  =  aK  F  F  F  (kjc)  •  Q(kJc)  •  F.  (-k.-k) 

m  aDL  L  mxy  x  y  me  x  y 


with  F^  given  by  Eq.  (59),  and  and  given  by  Eqs.  (35)  and  (36).  The  impedance  matrix 
elements  are  computed  using  Eq.  (51); 


4  V  y  F  (k  Jk  )  •  Q(k  Jc )  •  F  (-k  .  k  ) 

aO  L»  u  mxy  xy  nxy 

i  s  ~oo  j  a  -OO 


The  expansion  and  test  modes  F  and  F  are  exclusively  entire  domain  modes  given  by  Eqs. 

_  inc  n 

(55)  and  (56),  and  Q  is  given  in  the  Appendix.  It  was  found  that  using  the  r=(l,3,5) 
x-directed  modes  and  the  s=(2,4,6)  modes  on  each  patch  gives  convergent  results  for  E-plane 
scan.  For  off  E-plane  scans,  the  s=(4,6)  y-directed  modes  are  replaced  by  s=(  1 ,3)  to  account 
for  phase  asymmetry  in  the  y-direction. 

The  matrix  equation  is  then  solved  for  the  unknown  current  amplitudes,  and  the 
input  impedance  is  calculated  by  the  reaction  integral6; 


probe 


where  E  is  the  field  of  the  patch  currents,  and  the  impressed  probe  current  I  is 

scat  ‘  probe 

chosen  to  be  one  ampere.  Using  Eqs.  (7),  (34),  and  (41),  Eq.  (69)  can  be  rewritten  as; 


F7robe(k  Jc  ). 

inc  x  y 


(70) 


Convergent  results  were  obtained  by  using  101  terms,  (-50  <  i,  j  <  50),  in  each  of  the  series, 
and  evaluation  of  the  input  impedance  takes  approximately  3  minutes  on  a  DEC  MicroVax  III. 

With  the  input  impedance  known,  the  scanning  reflection  coefficient  can  be  computed 
by; 
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m<M  = 


(71) 


Z  (0,4»)  -  Z  (0  ,0  ) 

in _ i  n 

Z  (9,4>)  -  Z*  (0V) 

in  i  n 


when  me  array  is  conjugate  matched  at  broadside. 

The  idealized  probe  model  that  was  used  is  identical  to  the  one  used  in  an  earlier 
study.6  Although  the  model  ignores  the  self-reactance  of  the  probe  and  does  not  enforce 
continuity  of  current  at  the  probe-patch  intersection,  it  has  been  found  that  the 
approximation  works  well  for  substrates  that  are  less  than  0.05  dielectric  wavelengths  thick. 
For  larger  thicknesses,  the  self-reactance  of  the  probe  and  continuity  of  current  must  be 
accurately  modeled,  and  several  studies12,  3,14  have  addressed  these  details  carefully. 

Figure  5  shows  the  scanning  reflection  coefficient  for  a  type  A  element.  There  is  a 
total  reflection  at  0-57.9  degrees  in  the  E-plane,  and  no  blindness  in  the  H-plane.  As  with 
an  earlier  study6,  it  was  found  that  the  Green’s  function  for  y-directed  modes  contains  a 
surface  wave  resonance  in  the  H-plane  scan.  The  resonance  causes  the  Green’s  function  to  be 
singular,  and  a  steady  state  solution  does  not  exist  at  that  point.  The  value  of  the 
reflection  coefficient  at  that  angle  is  therefore  defined  to  be  the  value  at  a  vanishingly 
small  offset  angle,  R(0  +A0),  where  0  is  the  angle  of  the  surface  wave  resonance. 

8w  8W 

The  microstrip  line-coupled  elements  B  and  C  require  all  three  types  of  current 
expansion  modes.  The  electric  field  integral  equation  for  these  elements  can  be  written  as; 


where  R  is  the  unknow  reflection  coetficient,  s  is  the  surface  of  the  microstrip  line,  s  is 

M  n 

the  surface  of  the  n01  expansion  mode,  and  N  is  the  total  number  of  entire  domain  and 

piecewise  sinusoidal  modes.  The  first  integral  represents  the  fields  of  the  incident  and 
reflected  traveling  wave  modes  on  the  microstrip  feedline,  and  the  series  of  integrals 
represents  the  fields  of  the  entire  domain  and  piecewise  sinusoidal  modes. 

To  replicate  the  fields  of  a  traveling  wave  mode  on  the  microstrip  feed  line  in  an 


12.  Chew,  W.  and  Kong,  J  (1981)  Analysis  of  a  Circular  Microstrip  Disk  Antenna  with  a  Thick 
Dielectric  Substrate,  IEEE  Trans.  Antennas  and  Prop.,  AP-29:68-76. 

13.  Liu,  C.,  Hessel,  A.,  and  Shmoys.  J.  (1985)  Performance  of  Probe-Fed  Microslrip-Patch  Element 
Phased  Arrays,  Phased  Arrays  1985  Symposium  Digest,  pp.  157-176. 

14.  Aberle,  J.  and  Pozar,  D.  (1989)  Analysis  of  Infinite  Arrays  of  Probe-Fed  Rectangular 
Microstrip  Patches  Using  a  Rigorous  Feed  Model,  /EE  Proceedings,  136,  Pt.  H,  No.  2. 
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hi  ■  .159  cm 
h2  “  .159  cm 
XI  =  Y1  «  1.8  cm 
X2  =  Y2  =  1.95  cm 
C  ,  «  2.33-J.002 
Bg  =  2.33-J.002 

a  =  3.35  cm 
b  *  3.60  cm 


SCAN  ANGLE,  deg 


-0-  E-PLANE  -A-  H-PLANE 


Figure  5.  Computed  Reflection  Coefficient  Magnitude  of  Type  A  Element  as  a  Function 
of  Scan  Angle 
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infinite  array,  it  was  necessary  to  devise  a  method  of  extending  the  feed  lines  outside  of  the 
unit  cell  while  maintaining  the  infinite  array  periodicity.  On  physical  grounds,  one  can  argue 
that  the  traveling  wave  fields  are  bound  closely  to  the  microstrip  line,  and  will,  therefore, 
only  excite  currents  on  the  unit  cell  it  is  feeding.  This  approximation  will  be  good  for 
microstrip  lines  that  are  close  to  the  ground  plane,  and  which  do  not  come  in  close  proximity 
with  other  lines  or  patches.  Mathematically,  this  approximation  can  be  implemented  by  using 
the  single  element  Green’s  function  for  the  traveling  wave  field  integral  of  Eq.  (72),  and  the 
infinite  array  Green’s  function  for  the  series  of  integrals  representing  the  fields  due  to  the 
entire  domain  and  piecewise  sinusoidal  modes.  The  infinite  array  Green’s  function  is  used  for 
the  piecewise  sinusoidal  modes  because  they  account  for  the  non -propagating  modes  excited  on 
the  microstrip  line,  and  they  radiate  fields  that  are  strong  enough  to  excite  currents  on 
other  elements  of  the  array.  The  infinite  array  Green’s  function  must  also  be  used  for  the 
entire  domain  modes  because  they  represent  the  currents  on  the  microstrip  patches,  and  they 
couple  substantially  from  element  to  element.  Note  that  the  infinite  array  Green’s  function  of 
Eq.  (32)  reduces  to  the  single  element  Green’s  function  of  Eq.  (26)  if  there  is  no 
cell-to-cell  coupling. 

The  moment  method  procedure  described  earlier  is  then  applied  to  Eq.  (72)  to  form  a 
set  of  linearly  independent  equations.  An  additional  piecewise  sinusoidal  test  mode  provides 
the  extra  equation  needed  to  solve  for  the  N  +1  unknowns.  The  additional  unknown  comes 

ED+PWS 

from  the  unknown  reflected  traveling  wave  coefficient,  R.  The  matrix  equation  to  be  solved  for 
is  given  by; 


(73) 


where  the  subscript  R  refers  to  the  reflected  traveling  wave  expansion  mode. 

The  expressions  for  the  voltage  vector  and  impedance  matrix  elements  are  given  by; 


V 

m 


oo 


F  (-k  ,-k  )  •  Q(k  Jc )  ■ 

m  x  y  x  y 


F  (k  Jc  )  dk  dk 

i  nc  x  y  x  y 


.  oo  .  oo 


(74) 
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■ill 


F  (-k  ,-k  )  •  Q(k  Jc )  •  F™(k  Jc  )  dk  dk 

m  x  y  x  y  ref  x  y  x  y 


^mn  iH 

(  n*R) 


I  I 

|  =  -OO  j  =  .00 


F  (k  Jc  )  •  Q(k  Jc  )  •  F  (-k  ,-k  ) 

m  x  y  x  y  n'  x  y 


where  F.^  and  F^f  are  given  by  Eqs.  (61)  and  (62),  Q  is  given  in  the  Eqs.  (88)  -  (106) 

in  the  Appendix,  and  F^  and  F  are  the  piecewise  sinusoidal  and  entire  domain  expansion  and 
test  modes  given  by  Eqs.  (55),  (56),  and  (64).  As  with  the  probe  fed  patch,  it  was  found  that 
six  entire  domain  modes  on  each  patch  gave  convergent  results.  To  improve  the  modeling  of  'the 
currents  on  the  microstrip  line  in  the  vicinity  of  the  patch,  a  set  of  eight  piecewise 
sinusoidal  modes  are  used  to  compensate  for  the  non-propagating  currents  near  the 
discontinuities  for  elements  B  and  C.  For  element  B,  the  microstrip  traveling  wave  modes  and 
the  piecewise  sinusoidal  modes  overlap  the  patch  ,  and  it  was  found  that  the  results  vary  only 
slightly  for  overlaps  between  0.4  and  0.6  times  the  x-dimension  of  the  patch. 

The  numerical  integrations  in  Eqs.  (74)  and  (75)  are  carried  out  in  polar 
coordinates  (5  and  a,  defined  as; 

P2  =  k?  +  k?  (77) 


a  =  tan 


r  k  1 

hr)- 


The  limits  of  integration  are  P=0  to  °o  and  a=0  to  2k.  In  order  to  avoid  the  singularities  of 
the  spectral  Green’s  function  at  the  surface  wave  poles  and  the  traveling  wave  poles  on  the 
real  axis,  a  path  deformation  technique  is  used.15  By  introducing  an  imaginary  component  of 
p,  the  integration  path  can  be  deformed  sufficiently  far  away  from  the  real  axis  to  avoid 
numerical  difficulties,  as  shown  in  Figure  6.  The  modified  integration  path  will  give  the  same 
result  as  the  real  axis  integration  path  according  to  Cauchy’s  path  independence  rule,  which 


15.  Newman,  E.  and  Forrai,  D.  (1987)  Scattering  from  a  Microslrip  Patch,  IEEE  Trans.  Antennas 
and  Prop.,  AP-35. 
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I  MAG 


Figure  6.  Complex  Integration  Path  Deformation 


states  that  two  integration  paths  in  the  complex  plane  give  identical  results  as  long  as  no 

poles  are  enclosed  by  the  area  between  them.  Values  of  A=0.1k  ,  B=0.1k  ,  and  C=5.0k  are  used, 

0  0  0 

and  branch  cuts  are  chosen  such  that  the  values  of  k  represent  outward  traveling  and  decaying 
waves. 

The  integration  technique  was  checked  by  comparing  self-and-mutual  impedance  matrix 
elements  with  those  from  earlier  studies  on  printed  dipoles16,17,  and  the  values 
agreed  very  closely. 

With  the  reflection  coefficient  known  as  a  function  of  scan  angle,  the  scanning 
impedance  can  be  calculated  as; 

z»<e-«  ■ z.  !  -  •  <79> 


16.  Rana,  I.  and  Alexopoulos,  N.  (1981)  Current  Distribution  and  Input  Impedance  of  Printed 
Dipoles,  IEEE  Trans.  Antennas  and  Prop.,  AP-29:99-105. 

17.  Alexopoulos,  N.  and  Rana,  I.  (1981)  Mutual  Impedance  Computation  Between  Printed  Dipoles, 
IEEE  Trans.  Antennas  and  Prop.,  AP-29(No.  1):  106-1 11. 
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Using  a  DEC  MicroVax  III,  computation  of  the  input  impedance  for  a  type  B  element  takes 
approximately  ten  minutes,  and  for  type  C  elements  it  takes  approximately  eight  minutes. 

Figure  7  shows  the  scanning  reflection  coefficient  for  a  type  B  element.  There  is 
a  reflection  coefficient  magnitude  of  0.96  at  0=57.9  degrees  in  the  E-plane,  and  there  is  no 
blind  spot  in  the  H-plane.  Figure  8  shows  the  scanning  reflection  coefficient  of  a  type  C 
element.  There  is  a  reflection  coefficient  magnitude  of  0.8  at  47.3  degrees  in  the  E-planc, 
and  no  blindness  in  the  H-plane. 


2.6.  Mutual  Coupling  Analysis 


With  the  reflection  coefficient  known  as  a  function  of  scan  angle,  it  is  possible 
to  compute  the  mutual  coupling  coefficients  between  any  pair  of  elements  in  the  infinite 
array.  This  is  done  by  expressing  the  reflection  coefficient  as  a  Fourier  series  expansion; 

oo  oo 

R(u,v)  =  y  V  C  .-■imkau  e-j"kbv  (80) 

u  u  mn,00 

m  =  n  =  00 

where 


mn.OO 


k2ab 

4tc2 


R(u,v)  eik<ma“+nbv)  du  dv 


-TC  -71 
k  1  kb 


(81) 


is  the  expression  for  the  mutual  coupling  coefficient  between  the  center  element  and  the  m,nth 
element,  where  the  element  spacings  a  and  b  are  less  than  or  equal  to  one  half  of  a 
wavelength.  By  symmetry,  Eq.  (81)  can  be  reduced  to; 


mn ,  00 


7t  7C 
,k*  kb 

R(u,v)cos(mkau)cos(nkbv)  du. 

o  o 


k2ab 


Thus,  it  can  be  seen  that  the  mutual  coupling  coefficients  of  an  infinite  array  can  be 
obtained  from  the  scanning  reflection  coefficient,  R(u,v).  As  will  be  seen  in  the  next 


(82) 


18.  Hansen,  R.,  (ed.)  (1966)  Microwave  Scanning  Antennas,  Vol.II.  Array  Theory  and  Practice , 
Academic  Press,  New  York. 
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Figure  7.  Computed  Reflection  Coefficient  Magnitude  of  Type  B  Element  as  a  Function 
of  Scan  Angle 
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Figure  8.  Computed  Reflection  Coefficient  Magnitude  of  Type  C  Element  as  a  Function 
of  Scan  Angle 
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section,  this  method  can  be  used  to  approximate  the  scanning  reflection  coefficient  of  a 
finite  array  by  measuring  the  CmnQ0  coefficients,  and  computing  R(u,v)  from  Eq.  (80) 
with  a  truncated  series. 

2.7.  Scan  Blindness  Prediction 

Scan  blindness  is  a  condition  where  little  or  no  power  can  be  transmitted  or 
received  by  a  phased  array  for  certain  scan  angles,  as  indicated  by  the  large  reflection 
coefficient  magnitudes  in  Figures  5,  7,  and  8.  Equation  (80)  shows  that  the  reflection 
coefficient  is  due  to  the  mismatch  of  the  isolated  element  plus  contributic  as  from  all  other 
elements;  it  is  possible  for  an  in-phase  accumulation  of  this  coupled  power  to  cause  a  large 
reflection  at  one  or  more  scan  angles.  This  phenomenon  has  been  studied  for  several  infinite 
microstrip  array  geometries.5,6’7 

When  the  propagation  constant  of  a  surface  wave  in  the  loaded  dielectric  slab 
matches  a  Roquet  mode  propagation  constant,  it  has  been  found  that  a  forced  surface  wave  will 
be  excited,  such  that  no  real  power  enters  or  leaves  the  array  surface.  Mathematically,  this 
condition  can  be  expressed  as5; 

(p.„f  --  p?2-  +  ko“f  -  pir-  +  vf-  (83> 

It  has  also  been  found5  that  can  be  approximated  by  the  surface  wave  propagation  constant 
of  the  unloaded  dielectric  slab  over  a  ground  plane.  For  the  two-layer  dielectric  slab,  the  TM 
surface  wave  propagation  constants  correspond  to  zeros  of  the  transcendental  equation; 

— 'cos  k'h  fe  k  'cos  k  'h  +jk  'sin  k  'h  ] 
e(  I  2  2(_  1  o  1  1  J  1  1  ij 

+  jk  'sin  k  'h  [k  'cos  k  'h  +je  k  'sin  k  'h  1  =0  (84) 

J  2  2  2(  1  I  t  J  2  0  1  lJ 

where  k',  k-,  and  k'  are  defined  by  Eqs.  (19),  (20),  and  (21).  The  TE  surface  wave  propagation 
constants  are  given  by  the  zeros  of  the  equation; 

k  'cos  k  'h  (k  'cos  k  'll  +jk  'sin  k  'h  1  + 

2  2  2^1  I  1  J  0  1  lJ 

jk  'sin  k  Ti  fk  'cos  k  'h  +jk  'sin  k  'h  1  -  0.  (85) 

J  1  2  2(  0  1  1  J  1  1  'J 

It  can  be  seen  that  these  expressions  appear  in  the  denominator  of  the  two-layer  grounded 
dielectric  slab  Green’s  functions  given  in  the  Appendix. 
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An  earlier  study  points  out  that  leaky  wave  blindnesses  are  also  possible  in 
printed  circuit  arrays,  and  are  particularly  noticeable  when  the  patches  are  operating  at 
frequencies  away  from  resonance.  The  leaky  wave  blindness  angles  can  be  predicted  using  an 
approximation  similar  to  Eq.  (83)  with  replaced  by  P^,  the  propagation  constant  of  a 
leaky  wave  in  the  unloaded  dielectric  slab.  The  values  of  p(^  can  be  computed  from  the  complex 
zeros  of  Eqs.  (84)  and  (85),  with  both  real  and  imaginary  parts  of  the  solution  greater  than 
zero.  Note  that  the  Green’s  function  used  in  the  full  analysis  of  the  previous  sections  will 
account  for  all  surface  waves  and  leaky  waves  in  the  dielectric. 

The  blindspot  angles  predicted  by  Eq.  (83)  were  within  one-tenth  of  a  degree  of 
those  computed  by  the  full  wave  analysis  in  Figures  5,  7,  and  8. 

3.  EXPERIMENTAL  RESULTS 
3.1.  Waveguide  Simulator 

To  validate  the  theory,  a  rectangular  waveguide  simulator  was  constructed  using  the 
S-band  waveguide  shown  in  Figure  9.  The  lowest  order  TEjq  mode  in  the  waveguide  replicates  the 
fields  seen  by  an  element  in  an  infinite  periodic  array  radiating  a  plane  wave  in  the  H-pIane 
at  an  angle  ; 


e  =  Sin-'(-V]  (86) 

where  2b  is  the  major  dimension  of  the  waveguide.  The  waveguide  is  terminated  in  a  matched 
load  to  prevent  reflections  back  towards  the  radiating  element. 

A  basic  requirement  of  the  waveguide  simulator  is  that  the  element  must  have 
symmetry  about  its  center  on  the  x  and  y  axes.  Each  of  the  three  element  types  tested  have 
some  form  of  asymmetry.  As  can  be  seen  in  Figure  1,  element  A  has  a  probe  that  is  not 
symmetrical  about  the  center  of  the  x-axis,  and  elements  B  and  C  have  microstrip  feed  lines  on 
one  side,  also  violating  the  x-axis  symmetry.  However,  these  features  radiate  very  little  in 
comparison  with  the  microstrip  patches  and,  therefore,  will  not  substantially  excite  any 
waveguide  modes  that  would  cause  errors  in  the  simulation. 

All  simulator  measurements  were  made  on  a  calibrated  Hewlett  Packard  8510  network 
analyzer.  Element  A  was  fed  by  an  SMA  coaxial  probe  connector,  with  the  phase  reference  at  the 
ground  plane.  The  microstrip  lines  of  elements  B  and  C  were  fed  by  a  coax-to-microstrip 
transition  through  the  waveguide  wall,  with  the  phase  reference  at  the  transition. 

Figures  10-12  show  the  simulator  measurements  and  computed  data  for  the  three 
infinite  array  element  types.  Excellent  agreement  was  obtained  between  theory  and  experiment 
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Figure  9.  Waveguide  Simulator 
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Figure  10.  Measured  and  Computed  Input  Impedance  of  Type  A  Element  in  Waveguide  Simulator 
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Figure  11.  Measured  and  Computed  Input  Impedance  of  Type  B  Element  in  Waveguide  Simulator 
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Figure  12.  Measured  and  Computed  Input  Impedance  of  Type  C  Element  in  Waveguide  Simulator 
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in  all  three  cases.  Notice  the  looped  impedance  curves  of  elements  A  and  B  due  to  the  double 
resource...  The  second  small  loop  in  the  impedance  curve  of  Figure  10  between  4.8  and  5.0  GHz 
is  believed  to  be  due  to  a  small  amount  of  leakage  from  small  gaps  between  the  ground  plane  of 
the  elements  and  the  walls  of  the  waveguide.  Ideally,  those  edges  should  be  soldered  together, 
but  for  practical  purposes,  the  electrical  contact  between  the  element  ground  plane  and  the 
waveguide  walls  was  made  by  concocting  tape. 

The  impedance  characteristics  of  A  and  B  are  very  similar,  and  this  can  be  attributed 
to  the  fact  that  a  microstrip  feed  line  behaves  very  much  like  a  vertical  probe  feeding  the 
edge  of  a  microstrip  patch.19  Note  that  there  is  a  phase  rotation  between  the  impedance 
curves  of  elements  A  and  B  due  to  the  different  phase  references. 

Element  C  has  only  a  single  resonance  since  it  has  only  a  single  patch  and, 
therefore,  has  an  inherently  smaller  bandwidth.  As  can  be  seen  in  Figure  12,  the  predicted  and 
measured  results  are  nearly  identical.  'Hie  close  agreement  is  due  to  the  fact  that  the  type  C 
element  does  not  have  any  complicated  feed  current  behavior  to  model,  and  the  single  resonance 
curve  is  smoothly  varying. 

As  can  be  seen  from  Eq.  (86),  the  waveguide  simulator  cannot  replicate  the  Fields  of  an 
infinite  array  at  broadside  because  it  would  require  the  wavelength  to  go  to  zero.  In 
practice,  the  simulated  incidence  angle  varies  between  20  and  30  degrees,  depending  upon  the 
frequency. 

3.2.  11  X  11  Element  Array 

Although  the  waveguide  simulator  is  a  very  convenient  method  of  simulating  an  infinite  array, 

the  range  of  scan  angles  and  element  spacings  that  can  be  simulated  are  very  limited.  As  a 

separate  check  of  the  analysis,  an  1 1  x  1 1  array  of  type  A  elements  was  constructed,  and  the 

disassembled  array  is  shown  in  Figure  13.  It  has  been  found  that  five  match-loaded  elements  on 

all  sides  of  a  center  element  give  a  reasonably  close  approximation  to  an  piemen t  in  an 

infinite  array.  More  elements  will  be  required  to  simulate  the  occurrence  of  any  surface  wave 
20 

resonances. 

The  coupling  coefficients  between  the  center  element  and  all  other  elements  were 
measured  with  the  array  match-loaded.  Using  the  method  described  earlier,  the  mutual  coupling 
coefficients  were  computed  using  Eqs.  (71)  and  (82).  Figure  14  shows  the  measured  and  computed 


19.  Deshpande,  M.  and  Bailey,  M.  (1982)  Input  Impedance  of  Microstrip  Antennas,  IEEE  Trans. 
Antennas  and  Prop.,  AP-30(No.4):645-650. 

20.  Amitay,  N..  Galindo,  V.,  and  Wu,  C.  (1972)  Theory  and  Analysis  of  Phased  Array  Antennas , 
Wiley  Interscience.  New  York,  p.  366. 
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Figure  13.  11  x  11  Array  of  Type  A  Elements 
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values  for  one  quadrant  of  the  array  at  3.53  GHz.  The  results  agree  quite  well,  and  most  of 
the  values  are  closer  than  ldB  in  amplitude  and  10  degrees  in  phase.  The  larger  errors  occur 
for  values  of  coupling  less  than  35  dB,  as  would  be  expected. 

Using  the  measured  mutual  coupling  coefficients,  the  active  reflection  coefficient 
of  the  11  x  11  array  was  computed  by  Eq.  (80)  for  E-plane  and  H-plane  scans  from  0  degrees  to 
50  degrees.  Figure  15  shows  the  scanning  impedance  of  the  11  x  11  array  compared  to  the 
computed  infinite  array  scanning  impedance.  The  measured  and  computed  values  agree  quite  well 
at  the  lowest  frequency,  but  begin  to  depart  at  the  higher  frequencies.  There  appears  to  be  a 
disagreement  between  the  measured  data  and  the  computed  data  that  increases  with  frequency. 

The  most  likely  explanation  for  the  difference  is  that  the  finite  array  is  not  sufficiently 
large  enough  to  model  the  infinite  array.  As  mentioned  before,  a  previous  study  has  shown 
that  the  accurate  simulation  of  infinite  array  with  surface  wave  effects  requires  an  array 
larger  than  11  x  11. 

4.  BANDWIDTH  LIMITATION  STUDY 

In  the  interest  of  defining  the  limitations  of  proximity-coupled  patch  arrays,  a 
bandwidth  study  was  done  for  the  type  A  element.  To  begin  with,  the  bandwidth  of  an  uncovered 
probe-fed  square  patch  array  with  one-half  wavelength  spacing  was  computed  as  a  function  of 
dielectric  thickness.  The  resonant  patch  dimension  was  adjusted  so  that  the  input  impedance 
was  purely  real  at  broadside,  and  the  bandwidth  for  a  2:1  VSWR  was  calculated.  The  resulting 
bandwidth  for  the  uncovered  patches  is  shown  by  the  lower  solid  line  in  Figure  16.  The  upper 
solid  line  gives  the  location  of  the  E-plane  scan  blindness  angle  as  a  function  of  substrate 
thickness  for  one-half  wavelength  array  spacing. 

Next,  a  second  layer  of  proximity-coupled  patches  with  the  same  dielectric  constant 
as  the  first  layer  was  added,  and  the  upper  patch  dimension  was  adjusted  for  a  purely  real 
input  impedance.  The  broadside  bandwidth  was  computed  as  the  upper  layer  thickness  was 
increased.  The  dashed  curves  in  Figure  16  show  the  results  for  a  set  of  lower  layer 
thicknesses  ranging  from  0.01^o  to  0.05XQ.  As  the  upper  layer  thicknesses  are  increased,  it 
appears  that  the  bandwidths  increase  asymptotically  towards  a  limit  which  is  approximately  2.5 
times  the  uncovered  element  bandwidth. 

Paschen19  has  tabulated  the  bandwidth  improvement  factors  obtained  by  adding 
lossless  resonant  sections  to  a  complex  load.  The  values  were  calculated  from  the  Fano 

a 

gain-bandwidth  integral  constraint  given  by  ; 
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Figure  15.  Measured  and  Computed  Scanning  Impedance  of  11  x  11  Array 
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In  |  Rn(w)  |  do)  <  K  (87) 

o 

where  R^fio)  is  the  reflection  coefficient  of  the  system  with  n  resonant  matching  sections,  io 
is  the  radian  frequency,  and  K  is  a  constant  determined  by  the  impedance  of  the  complex  load. 
This  equation  implies  that  there  is  a  physical  limit  to  the  frequency  bandwidth  over  which  the 
reflection  coefficient  of  the  load  can  be  kept  below  a  certain  value.  The  bandwidth 
improvement  factors  for  n  resonant  sections  were  computed  to  be21; 

21.  Paschen,  D.  (1983)  Broadband  Microstrip  Matching  Techniques,  Allerton  Antenna  Applications 
Symposium  Digest ,  pp.1-20. 
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number  of  sections  bandwidth  improvement  factor 

1  2.33 

2  2.84 

3  3.18 

«  3.86 

The  bandwidth  improvement  of  2.33  for  a  single  resonant  matching  section  is  close  to 
the  bandwidth  improvement  of  2.5  obtained  by  adding  the  second  layer  of  resonant  patches  to 
the  uncovered  patch  array.  If  the  upper  patch  behaves  like  a  resonant  circuit  in  parallel  with 
the  lower  patch,  then  the  bandwidth  improvement  should  be  governed  by  the  Fano  gain-bandwidth 
integral  constraint.  Since  the  upper  patch  radiates,  it  will  look  like  a  lossy  resonant 
section,  and  will,  therefore,  make  the  bandwidth  improvement  slightly  larger  than  the  amount 
predicted  by  Fano  for  lossless  matching  sections.  From  a  practical  point  of  view,  this  means 
that  the  bandwidth  is  constrained  by  the  uncovered  lower  patch  array. 

As  a  check  of  these  predictions,  a  set  of  six  measurements  were  made  using  the 
infinite  array  simulator.  Figure  17  shows  the  different  layer  combinations  for  each  of  the  six 
cases.  The  upper  layer  thicknesses  were  chosen  to  study  the  effect  on  bandwidth  of  increasing 
thickness.  For  case  1,  the  lower  layer  thickness  was  fixed  at  0.159  cm.  Figure  18  shows  the 
measured  and  computed  input  impedance  for  case  1A,  which  is  an  uncovered  microstrip  element  in 
an  infinite  array.  The  bandwidths  were  computed  by  normalizing  the  impedance  to  the  purely 
real  input  impedance  where  the  curve  crosses  the  real  axis  of  the  Smith  chart.  Note  that  the 
scan  angle  of  the  array  shifts  a  small  amount  with  frequency  as  shown  by  Eq.  (86).  Therefore, 
the  bandwidth  is  computed  and  measured  as  the  array  scans  with  frequency  from  31  degrees  at 
4.0  GHz  to  20  degrees  at  6.0  GHz.  Figure  19  shows  the  measured  and  computed  input  impedance 
for  case  IB,  which  has  an  upper  layer  thickness  of  0.159  cm.  There  is  excellent  agreement,  and 
the  bandwidth  has  increased  by  a  factor  of  2.2  from  case  1A.  Figure  20  shows  the  measured  and 
computed  input  impedance  for  case  1C,  which  has  an  upper  patch  layer  thickness  of  0.318  cm. 

The  results  agree  very  closely,  and  the  bandwidth  has  increased  to  3.0  times  the  uncovered 
array  bandwidth. 

For  case  2,  the  lower  layer  thickness  was  fixed  at  0.318  cm.  Figure  21  shows  the 
measured  and  computed  input  impedance  of  case  2A,  which  is  an  uncovered  array.  Figure  22  shows 
the  measured  and  computed  input  impedance  for  case  2B,  which  has  an  upper  patch  layer 
thickness  of  0.159  cm.  The  bandwidth  has  increased  considerably  to  2.7  times  the  uncovered 
array  bandwidth.  Figure  23  shows  the  measured  and  computed  input  impedance  for  case  2C,  which 
has  an  upper  patch  layer  thickness  of  0.138  cm.  Although  the  upper  layer  thickness  has  been 
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—  upper  patch  —  lower  patch 


CASE-H 


Figure  17.  Substrate  Combinations  for  Bandwidth  Limitation  Study 
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XI 


hi  =  Ocm 
h2  =  .159  cm 
XI  =  Y1  =0cm 
X2  =  Y2=  1.95  cm 
=  2.33  -  j.002 
e2  =  2.33  -  j.002 


Measured  • 
(2:1  B.W.  3.9%) 


Figure  18.  Measured  and  Computed  Input  Impedance  for  Case  1A  of  Bandwidth  Limitation  Study 


X2 


Predicted  a 
(2:1  B.W.  3.7%) 
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CASE  IB: 

hi  =  .159  cm 
h2  =  .159  cm 
XI  =  Y1  =  1.8  cm 
X2  =  Y2  =  1.95  cm 
e1  =  2.33  -j. 002 
e2  =  2.33  -  j.002 


Measured  ®  Predicted  A 

(2:1  B.W.  8.5%)  (2:1  B.W.  9.0%) 


Figure  19.  Measured  and  Computed  Input  Impedance  for  Class  IB  of  Bandwidth  Limitation  Study 


CASE  1C: 


hi  =  .318  cm 
h2  =  .159  cm 
XI  =  Y1  =  1.7  cm 
X2  =  Y2  =  1 .95  cm 
6!  =  2.33  -j. 002 
e2  =  2.33  -  j.002 


XI 

I  h — hh 

X2 


i 

T 


hi 


Measured  ®  Predicted  ^ 

(2:1  B.W.  11.9%)  (2:1  B.W.  11.4%) 


Figure  20.  Measured  and  Computed  Input  Impedance  for  Case  1C  of  Bandwidth  Limitation  Study 
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CASE  2C: 


hi  =  .318  cm 
h2  =  .318  cm 
XI  =  Y1  =  1.7  cm 
X2  =  Y2  =  1 .95  cm 
Ei  =  2.33  -  j.002 
e2  =  2.33  -  j.002 


Measured  •  Predicted  A 

(2:1  B.W.  19.5%)  (2:1  B.W.  18.3%) 


Figure  23.  Measured  and  Computed  Input  Impedance  for  Case  2C  of  Bandwidth  Limitation  Study 
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doubled,  the  bandwidth  has  only  improved  a  small  amount  to  2.8  times  the  uncovered  array 
bandwidth.  The  results  can  be  summarized  as; 


case 

upper  thickness 

bandwidth 

bandwidth  improvement 

1A 

o.ox 

o 

3.9% 

1.0 

IB 

0.025A, 

o 

8.5% 

2.2 

1C 

0.05X 

o 

11.9% 

3.0 

2A 

o.ox 

0 

6.9% 

1.0 

2B 

0.025X 

0 

18.7% 

2.7 

2C 

0.05A, 

0 

19.5% 

2.8 

These  results  show  that,  for  these  parameters,  the  bandwidth  improvement  does  approach  a 
limit  of  2.5  to  3.0  times  the  uncovered  array  bandwidth  as  the  upper  patch  layer  thickness  is 
increased.  Note  that  these  results  represent  a  trend  for  the  particular  dielectric  constants 
and  layer  thicknesses  that  were  studied.  Changes  in  these  parameters  may  yield  better  or  worse 
bandwidth  improvements,  and  any  variations  must  be  analyzed  separately. 

5.  ARRAY  ELEMENT  DESIGN  STUDY 

The  bandwidths  of  the  antennas  described  in  the  previous  section  were  computed  at 
broadside,  and  will  be  reduced  by  the  input  impedance  variation  as  the  array  is  scanned  off  of 
broadside.  To  investigate  the  effects  of  scanning  upon  bandwidth,  the  scanning  impedance  of 
several  array  elements  were  analyzed. 

To  begin  with,  an  array  element  was  designed  for  a  broadside  bandwidth  of  18 
percent.  Figure  24  shows  the  input  impedance  of  the  element  normalized  to  50Q.  The  small 
circle  at  the  center  of  the  plot  represents  the  2:1  VSWR  region.  The  array  impedance  stays 
within  the  2:1  VSWR  circle  for  scans  out  to  60  degrees  in  the  E-plane  and  H-plane  at  f=0.91f0 
and  f=1.05fo,  but  a  serious  mismatch  occurs  at  the  center  frequency  f=  1  -Of0  in  the  H-plane. 

This  example  shows  that  the  scanning  impedance  variation  can  have  a  serious  effect  upon  the 
operating  bandwidth  of  an  array,  and  must  be  accounted  for  in  the  design. 

A  second  array  element  was  designed  for  scanning  bandwidth,  and  the  resulting  scan 
performance  is  shown  in  Figure  25.  The  array  has  a  14  percent  bandwidth  at  broadside,  and  can 
scan  over  a  50  degree  cone  with  a  10  percent  bandwidth.  This  result  was  achieved  by  reducing 
the  H-plane  spacing  of  the  elements  from  0.5X  to  0.45\Q,  and  the  aspect  ratio  of  the  patches 
from  1.2  to  1.0.  These  changes  reduced  the  H-plane  scanning  impedance  variation  enough  to 
maintain  a  match  at  the  center  of  the  band  at  the  expense  of  overall  bandwidth.  It  has  been 
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Figure  25.  Scanning  Impedance  of  Array  Element  Optimized  for  Scanning  Bandwidth 
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22 

suggested  that  a  low  £f  for  the  upper  layer  would  further  improve  the  scanning  bandwidth. 

6.  CONCLUSIONS 

A  full  wave  analysis  was  applied  to  three  types  of  proximity-coupled  microstrip  array 
elements.  The  infinite  array  spectral  dyadic  Green’s  function  was  used  so  that  all  mutual 
coupling  and  surface  wave  effects  were  included.  A  method-of-moments  was  used  to  solve  the 
electric  field  integral  equation,  and  three  types  of  expansion  modes  were  used  to  efficiently 
model  the  current  distribution  on  the  feed  line  and  patches.  By  using  the  approximation  that 
the  microstrip  feed  lines  couple  only  to  the  unit  cell  they  are  feeding,  the  incident  and 
reflected  quasi-TEM  fields  can  be  found  by  integrating  over  the  isolated  semi-infinite  line 
with  the  single  element  Green’s  function.  This  makes  it  possible  to  model  the  field  of  a 
traveling  wave  mode  with  a  source  current  that  extends  outside  of  the  unit  cell  while 
maintaining  the  infinite  array  periodicity. 

The  scan  performance  of  three  proximity-coupled  element  types  was  studied.  It  was 
found  that  scan  blindness  does  occur  due  to  surface  wave  resonances,  and  a  complete  E-plane 
blindness  was  found  for  element  A,  and  a  large,  non-unity  reflection  was  found  for  elements  B 
and  C. 

A  bandwidth  limitation  study  showed  that  the  bandwidth  improvement  factor  obtained 
by  adding  a  layer  of  proximity  coupled  elements  to  an  uncovered  patch  array  was  approximately 
2.5.  This  limit  may  be  related  to  the  theoretical  limit  for  broadband  matching  of  a  complex 
load  impedance,  as  determined  by  the  Fano  gain-bandwidth  integral  constraint.  Experimental 
measurements  agree  with  the  results  of  the  bandwidth  limitation  study. 

An  array  element  design  study  showed  that  the  input  impedance  of  a  proximity- 
coupled  array  element  can  vary  substantially  as  a  function  of  scan  angle.  It  was  found  that 
the  scan  performance  and  bandwidth  can  be  traded  off,  and  an  array  element  was  designed  with  a 
10  percent  bandwidth  over  a  50  degree  scan  cone.  The  results  of  this  study  imply  that  the 
design  of  a  proximity-coupled  microstrip  array  with  a  wide  scanning  bandwidth  must  incorporate 
mutual  coupling  and  surface  wave  effects. 


22.  Hansen,  R.  (1988)  in  personal  communications  with  the  author. 
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Appendix:  Spectral  Dyadic  Green’s  Function 


The  spectral  dyadic  Green’s  functions  that  are  needed  for  the  analysis  are  given 
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where 


p  =  field  direction  (x,y,  or  z) 
q  =  source  direction  (x,y,  or  z) 
r  =  1  for  source  at  z=h  +h2 
r  =  2  for  source  at  z=h‘ 
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(105) 


k  '2=  e .  k*  -  p2 

2  2  0 

(106) 

(107) 
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0 

(108) 

k  _  2n 
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The  z-direcied  components  were  obtained  by  integrating  over  the  z-directed  probe  of  constant 
current  and  height  h2. 
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